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ABSTRACT 



In a wireless communication system, radio frequency signal 
propagation characteristics associated with a signal trans- 
mitted from several locations in the coverage area of the 
wireless communication system are measured. Next, a plu- 
rality of signature regions arc defined in the coverage area, 
wherein the signature regions each contain several selected 
locations having selected measured radio frequency signal 
propagation characteristics. For a selected signature region, 
the associated measured radio frequency signal propagation 
characteristics are processed to produce a spatially averaged 
radio frequency signature for the selected signature region. 
The spatially averaged radio frequency signature represents 
the dominant features of the radio frequency signal propa- 
gation characteristics of a signal transmitted within the 
selected signature region. Finally, an association between 
the spatially averaged radio frequency signature and the 
selected signature region is recorded. 

12 Claims, 8 Drawing Sheets 



o-( start) 



202Hsrlect a location I 



FOE SELECTED LOCATION. MEASURE 
BAD (0 FREQUBNCT SIGNAL PROPAGATION 
CHARACTERISTICS Of SIGNAL FROM TRANSMITTER 



1 



^NBIT LOCATION?^ — (SELECT NEXT LOCATION | 
206 \H0 



DEFINE A PLURALITY OP BIGNATUBE 
REGIONS, EACH CONTAINING A PLURAL ITT 
07 SELECTED LOCATIONS 



SELECT A SIGNATURE BECI0N I 



FOB A SELECTED SIGNATURE REGION, PROCESS THE MEASURED 
RADIO FREQUENCY SIGNAL PROPAGATION CHARACTERISTICS 
WITBIN TBI SELECTED SIGNATURE REGION TO PRODUCE A 
RADIO FREQUBNCT SIGNATURE FOR THB SELECTED SIGNATURE 
REGION THAT REPRESENTS TUB DOMINANT FEATURES OF THE 
RADIO FREQUENCY SIGNAL PROPAGATION CHARACTERISTICS 
IN TBE SELECTED SIGNATURE REGION 



RECORD AN ASSOCIATION BETWEEN THE 
RADIO FREQUENCY SIGNATURE AND THE 
SELECTED SIGNATURE REGION 



<NEIT SIGNATURE BEGlONt^ r—* 
218 ~W 
222-C END } 



SELECT NEXT 
SIGNATURE REGION 

1 " 

220 



03/10/2004, EAST version: 1.4.1 




03/10/2004, EAST Version: 1.4.1 



U.S. Patent Jul. 3, 2001 Sheet 2 of 8 US 6,256,506 Bl 



200 - (START) 
202 H SELECT A LOCATION] 



FOR SELECTED LOCATION, MEASURE 
204 -H RADIO FREQUENCY SIGNAL PROPAGATION 

CHARACTERISTICS OF SIGNAL FROM TRANSMITTER 



<NEXT LOCATION?> — H SELECT NEXT LOCATION | 



206 



210 ~ 



208 



NO 



DEFINE A PLURALITY OF SIGNATURE 
REGIONS, EACH CONTAINING A PLURALITY 
OF SELECTED LOCATIONS 



212' SELECT A SIGNATURE REGION 



FOR A SELECTED SIGNATURE REGION, PROCESS THE MEASURED 
RADIO FREQUENCY SIGNAL PROPAGATION CHARACTERISTICS 
WITHIN THE SELECTED SIGNATURE REGION TO PRODUCE A 
RADIO FREQUENCY SIGNATURE FOR THE SELECTED SIGNATURE 
REGION THAT REPRESENTS THE DOMINANT FEATURES OF THE 
RADIO FREQUENCY SIGNAL PROPAGATION CHARACTERISTICS 
IN THE SELECTED SIGNATURE REGION 



RECORD AN ASSOCIATION BETWEEN THE 
216 H RADIO FREQUENCY SIGNATURE AND THE 
SELECTED SIGNATURE REGION 



^NEXT SIGNATURE REGION?^ 



YES 



218 



NO 



-214 



222^ END ) 



SELECT NEXT 
SIGNATURE REGION 

\ 

220 



FIG. 2 



03/10/2004, EAST version: 1.4.1 



U.S. Patent Jul. 3, 2001 Sheet 3 of 8 US 6,256,506 Bl 



s(i+3) 



300 



304 



xxxxxxxxxxxxxxxxxxxxxxxxxxxx 



302 r 



TIME 



S(i)= 



V i+N-1 



TIME 



F/G. 3 



s(i) 



IMPULSE RESPONSE ESTIMATOR |— » hfi> 



J06 



h 



G)m J 



TIME 
DELAY 

i 



500- ( start) 



502- 



SELECT FIRST TIME | 



504- 



RECALL MEASURED CHANNEL IMPULSE 
RESPONSE(S) FOR SELECTED TIME 



506- 



VECTORIZE RECALLED IMPULSE RESPONSE(S) 
AND APPEND AS COLUMN IN MATRIX G, WHEREIN 
COLUMN CORRESPONDS TO SELECTED TIME 



508* 



^ALL TIMEs"?y ^-» [sE~LECT NEXT TIME | -l 



5f2- 



510 



YES 



EXTRACT PRINCIPLE SUBSPACE, S, OF MATRIX G 

R=GG H 

R=UAU H 
U=[ui u 2 »-u N ] 

A«diag([X! X 2 — X N J) 
Up=[ui U2-"Up], P£N 
Ap 8 diag([Xi X2— Xp]) 
S=UpA P 



514 -( END ) 



F'lG-. 5 



03/10/2004, EAST version: 1.4.1 



U.S. Patent 



Jul. 3, 2001 



Sheet 4 of 8 



US 6,256,506 Bl 



400- ( START) 
402 -" j RECEIVE LOCATION REQUEST | 



I 



MEASURE RADIO FREQUENCY SIGNAL 
404-1 PROPAGATION CHARACTERISTICS 
OF SIGNAL FROM TRANSMITTER 



I 



/ K YES I 1 

^NEED ANOTHER MEASUREMENT? > »| DELAY [^408 



406 



NO 



PROCESS THE MEASURED RADIO FREQUENCY SIGNAL 
PROPAGATION CHARACTERISTICS TO PRODUCE A RADIO 
FREQUENCY SIGNATURE FOR THE SIGNAL FROM THE 
SUBSCRIBER UNIT TO BE LOCATED THAT REPRESENTS THE 
DOMINANT FEATURES OF THE RADIO FREQUENCY SIGNAL 
PROPAGATION CHARACTERISTICS 



-4f0 



SELECT A DATA RECORD FROM THE SIGNATURE DATABASE 
HAVING DATA THAT REPRESENTS THE DOMINANT FEATURES 
OF THE RADIO FREQUENCY SIGNAL PROPAGATION 
CHARACTERISTICS OF A SELECTED SIGNATURE REGION 



-472 



COMPARE THE RADIO FREQUENCY SIGNATURE FOR THE 
SIGNAL FROM THE SUBSCRIBER UNIT TO BE LOCATED TO THE 
SELECTED DATA RECORD FROM THE SIGNATURE DATABASE 
TO PRODUCE A FIGURE OF MERIT 



416- 



■414 



STORE FIGURE OF MERIT FOR THE SELECTED DATA 
RECORD FROM THE SIGNATURE DATABASE 



^NEXT DATA RECORD?^ 
418 



YES 



422~ 



NO 



SELECT NEXT 
DATA RECORD 



RANK THE DATA RECORDS BASED ON THE ASSOCIATED 
FIGURES OF MERIT TO FIND A LOCATION ESTIMATE 



424- 



420 



OUTPUT A SIGNATURE REGION | 



426- ( END ) 



03/10/2004, EAST Version: 1.4.1 



U.S. Patent 



Jul. 3, 2001 



Sheet 5 of 8 



US 6,256,506 Bl 




03/10/2004, EAST Version: 1.4.1 



U.S. Patent Jul.3,2001 Sheet 6 of 8 US 6,256,506 Bl 



700H START 




702' SELECT FIRST SELECTED LOCATION 



704' 



1. 



RECALL MEASURED CHANNEL IMPULSE 
RESPONSE(S) FOR SELECTED LOCATION 



VECTORIZE RECALLED IMPULSE RESPONSE(S) 
706H AND APPEND AS COLUMN IN MATRIX G, WHEREIN 
COLUMN CORRESPONDS TO SELECTED LOCATION 



708 



I 



^ALL LOCATIONS?^ 



NO 



710 



SELECT NEXT LOCATION 



YES 



712- 



EXTRACT PRINCIPLE SUBSPACE, T, OF MATRIX F 

R»FF H 

R=UAU H 
U=[U! u 2 --u N ] 

A«diag([Xi X2—M) 
Up s [ui U2"«up], P£N 
Ap=diag([Xi X2 ••♦Xpl) 
T'UpAp 

I 



7M- ( END~) 



FIG. 7 



03/10/2004, EAST version: 1.4.1 



U.S. Patent Jul.3,2001 Sheet 7 of 8 US 6,256,506 Bl 




FIG. 9 



03/10/2004, EAST version: 1.4.1 



U.S. Patent Jul.3,2001 Sheet 8 of 8 US 6,256,506 Bl 




INDEX 




03/10/2004, EAST version: 1.4.1 



US 6,256 ; 

1 

METHOD AND SYSTEM FOR CREATING A 
RADIO FREQUENCY SIGNATURE FOR A 
SIGNATURE REGION IN A COVERAGE 
AREA OF A WIRELESS COMMUNICATION 

SYSTEM 5 

RELATED APPLICATIONS 

This application is related to application Ser. No. 09/475, 
095 entitled "METHOD AND SYSTEM FOR COMPAR- 
ING MEASURED RADIO FREQUENCY SIGNAL 
PROPAGATION CHARACTERISTICS IN A WIRELESS 
COMMUNICATION SYSTEM," which is filed on even 
date herewith, and which is incorporated herein by refer- 
ence. 35 

FIELD OF THE INVENTION 

The present invention is related in general to wireless 
communication systems, and more particularly to an 
improved method and system for creating a radio frequency 20 
(RF) signature for a signature region in a coverage area of 
a wireless communication system. 

BACKGROUND OF THE INVENTION 

For both users and operators of wireless communications 
systems it is often desirable to locate subscriber units that are 
making calls, or are otherwise transmitting or receiving data. 
Applications for technology that can locate subscribers 
include locating subscribers requesting 911-emergency 3Q 
services, assisting law enforcement agencies in various law 
enforcement activities, providing users with local area maps 
or directions, and providing users with advertising informa- 
tion related to the user's location. 

Several technologies have been either used or proposed 35 
for providing the subscriber unit's location in a communi- 
cation system, and each of these technologies has disadvan- 
tages. For example, a time difference of arrival (TDOA) 
system uses measurements from a plurality of base trans- 
ceivers at known locations. This type of system requires 4 q 
signal reception of the subscriber unit's signal at more than 
one base transceiver, time synchronization between the base 
transceivers, and the measurement of line-of-sight wave 
front propagation, without confusing the line-of-sight sig- 
nals with multipath signals, which are signals bounced and 45 
diffracted from objects affecting the propagation of the 
signal. The accuracy of synchronization of base transceivers 
in a cellular system based on interim specification 95 
(IS-95), a code division multiple access (CDMA) standard, 
is not accurate enough to make signal propagation measure- 50 
ments that will provide location estimates with the desired 
resolution. Additionally, in a dense downtown area, also 
known as the urban canyon, detection of the line-of-sight 
wave front may not be possible because line-of-sight signal 
components may be shadowed or obscured by close-in 55 
multipath components. 

More advanced TDOA systems use smart antennas to help 
eliminate some of the multipath components, enhancing 
detection of the prompt wave front. Smart antennas may also 
be able to add angle-of-arrival (AOA) measurements to the 60 
geometric equations for location finding, which further 
enhances performance. Unfortunately, the multipath envi- 
ronment of some urban canyons is too confused for smart 
antennas to be an effective solution. Smart antenna overlap 
solutions also require adding special phased array antennas 65 
to work in conjunction with the existing antennas. Service 
providers will most likely be reluctant to invest in additional 
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antennas until smart antennas are needed for capacity 
requirements as well as location finding requirements. 

Some engineers and operators propose using the global 
positioning system satellite system for determining a sub- 
scriber's location. This technology also has a problem in the 
urban canyon. In many places it is not possible to receive 
line-of-sight signals from a number of satellites that is 
needed to produce an accurate location estimate. 

Other proposals include placing many microcells in an 
urban core area so that the subscriber's location is estimated 
by identifying what cell site the subscriber is using. 
Although microcells are gaining popularity in certain areas, 
the market has generally rejected microcells for capacity 
coverage, which means that the number of microcells 
deployed is probably not extensive enough to support loca- 
tion finding. Most current communication deployment plans 
still specify macrocellular coverage. 

Therefore, it should be apparent that there remains a need 
for an improved, cost-efficient approach for locating sub- 
scriber units in a cluttered area of a communications system 
service area, where radio frequency signal characteristics of 
indirect or non -line-of-sight propagation signals are benefi- 
cially used to locate a subscriber unit. Accordingly, there is 
also a need for an improved method and system for creating 
a radio frequency signature region in a coverage area of a 
wireless communication system, wherein the radio fre- 
quency signature is responsive to the dominant features of 
radio frequency propagation characteristics of a signal trans- 
mitted from a selected signature region. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The novel features believed characteristic of the invention 
arc set forth, in the appended claims. The invention itself, 
however, as well as a preferred mode of use, further objects, 
and advantages thereof, will best be understood by reference 
to the following detailed description of an illustrative 
embodiment when read in conjunction with the accompa- 
nying drawings, wherein: 

FIG. 1 illustrates multipath signal propagation in an urban 
canyon area of a wireless communication system coverage 
area; 

FIG. 2 is a high-level logic flowchart that illustrates the 
method and operation of creating a spatially averaged radio 
frequency signatures for signature regions in a coverage area 
of a wireless communication system in accordance with the 
method and system of the present invention; 

FIG. 3 depicts the processing of raw measurement data to 
produce an implicit measurement of radio frequency propa- 
gation characteristics in the form of a channel impulse 
response measurement, which measurement is used to pro- 
duce database record representing a radio frequency signa- 
ture in accordance with the method and system of the present 
invention. 

FIG. 4 is a high-level logic flow chart that illustrates the 
method and operation of measuring a time averaged radio 
frequency signature for a subscriber unit in a coverage area 
of a wireless communication system in accordance with the 
method and system of the present invention; 

FIG. 5 is a more detailed logic flowchart of the process of 
producing a spatially averaged radio frequency signature, as 
shown in FIG. 4; 

FIG. 6 depicts a system used to create radio frequency 
signatures for signature regions in a coverage area of a 
wireless communication system, and also for measuring a 
radio frequency signature for a subscriber unit to be located 
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in a coverage area during a location event, in accordance metric decomposition of the received signal. However, in a 

with the method and system of the present invention; preferred embodiment, the radio frequency signal propaga- 

FIG. 7 is a more detailed logic flowchart of the process of tion characteristics are measured implicitly by sampling 

producing a spatially averaged radio frequency signature, as f ° ra P u " n S ^d recording the channel [impulse responses , of 

h FIG 2* jo ^ the particular channel between subscnber 38 and each of the 

s own in . , one Qr more antennas at Dase transceiver 40. 

FIG. 8 is a graphical representation of evolving power , f the fadio frequency signa i propagation characteristics 

delay profiles, which are used by the method and system of are measure d explicitly, a parametric model of the channel 

the present invention; impulse response is formed, as represented by the following 

FIG. 9 is a graphical representation of nonlinearly aver- 1Q equation: 
aging the data of FIG. 8; 

FIG. 10 is a graphical representation of power delay equation l 

profiles collected from a plurality of locations within a MO = 2j w « ( ' ,r ' ( ^ w '^ (,) *^ (,)>5( '" T « (/)) 
signature region; and 

FIG. 11 is a graphical representation of nonlinearly aver- is 

aging the data of FIG. 10. where: 

e is an index to the antenna of interest, 

DETAILED DESCRIPTION OF THE Q(l) repre sents the number of paths, 

INVENTION a ^ * s tne cortex amplitude of path q, 

With reference now to FIG. 1, there is depicted commu- 20 r r (8,*,Q is the complex valued response of antenna e due 

nications system coverage area 20, which includes buildings to stimulus by a wavefront with azimuth angle 6, 

22-30 that reflect and diffract radio frequency signals that elevation angle <fc and polarization angle |, and 

travel propagation paths 32-36 between subscriber unit 38 >s the lime delay of path q. 

and base transceiver 40. e 9 uat i°n for the nonstationary impulse response can 

' „ T ~ . . . it _ 25 be further decomposed by exploiting relative stationarity of 

It should be evident from FIG. 1 that propagation paths certam parts . 

32-34 have different times of flight and angles of arrival 

than propagation path 36. Thus, if transceiver 40 measures ^ equation 2 

the angle of arrival and timing information from propagating ^ r) = y a q {T)r e (0 q (n ?,(n fi,tT))tf(r - r,(D) 

signals and uses that data to estimate the location of sub- 30 pi 
scriber unit 38, the data from propagation paths 32-36 will 
not coalesce to indicate a single location in communication 

system coverage area 20 unless signal reflections and dif- In equation 2, T represents time and t represents time 
fractions caused by buildings 22-30 are taken into account. dela y- These parameters are sufficient to characterize the 
„ r . - . ™„ - 4 . . . . # « ... location dependent multipath distortion observed by the 
With reference now to HG. 2, there is depicted a high- 35 transmitted signal. At any given time, the transmitter will 
level logic flowchart that illustrates the method and opera- Qccupy a bcation> e g ( x (t),y(t)), for which the parameters 
tion of creating spatially averaged radio frequency signa- of the measurcd ^puise response are the radio frequency 
tures for signature regions in a coverage area in a wireless signal propa gation characteristics. (Note that this description 
communication system in accordance with the present 0 f the p arametr i c mo del excludes modeling of the Doppler 
invention. As illustrated, the process begins at block 200, 40 frequency shifts due to speed of travel through the environ- 
and thereafter passes to block 202 wherein the process ment. This data is implicitly contained within a^T). Dop- 
sclccts a location in communication system coverage area 20 pier shifts are due to Lravel speed, which is only secondarily 
(See FIG. 1). In a preferred embodiment, this process of associated with position. The other parameters described are 
selecting a location includes noting the present location of a more directly correlated with the transmitter's location, 
transmitter in communication system coverage area 20. This 45 However, the addition of Doppler could, in some 
may be accomplished by recording a location from a GPS circumstances, prove beneficial.) Implicit measurements of 
receiver located with the transmitter, or noting the location these characteristics are made by recording samples of the 
using maps or surveying equipment. Other electronic loca- impulse responses without reduction to the fundamental 
tion finding devices may also be used, such as the LORAN parameters. For example, a vector-valued implicit 
location finding system, operated by the United States Coast 50 measurement, h tf (T), of these radio frequency signal char- 
Guard. In a preferred embodiment, the selected location is acteristics can be represented as follows: 
one of a plurality of locations selected while driving or 

walking through communication system coverage area 20 AXO-IWO) hjr,T t ) . . . A e (r,(iV-i)r a )J r equation 3 

while carrying a transmitter that emits a signal. Some Note that any linear transformation of the responses (e.g. 

selected locations are shown in FIG. 1 at the symbol 55 Fourier transformation) yields an equally valid basis for 

indicated by reference numeral 42. implicit description of the radio frequency signal propaga- 

After a location is selected, the process measures radio tion characteristics. Certain transformations will create 

frequency signal propagation characteristics of the signal opportunities to limit the data to a smaller region of support, 

from the transmitter at the selected location, as indicated at thereby discarding unimportant data and reducing the 

block 204. The measured radio frequency signal propagation 60 amount of data to be recorded. For example, a Fourier 

characteristics may include the number of multipath transformation followed by band-limiting to the frequency 

components, the average power in each multipath bins in which there is significant signal power may reduce 

component, the gain and phase of the antenna due to the the amount of data recorded. These transformations follow 

azimuth and/or elevation angles of arrival and/or polariza- the methods of van der Veen, Vanderveen and Paulraj, "Joint 

tion of each component, the Doppler frequency of each 65 Angle and Delay Estimation Using Shift-Invariance 

component, and the relative time delay of each component. Techniques," IEEE Transactions on Signal Processing, Feb- 

These characteristics may be measured explicitly by para- ruary 1998. 
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In a preferred embodiment, RF signal propagation char- radio frequency signature represents the dominant features 

acteristics are measured implicitly by sampling voltage of the radio frequency signal propagation characteristics of 

measurements from one or more antennas at base transceiver a signal transmitted from the selected signature region and 

40. With reference now to FIG. 3, a single voltage sample is received at a base transceiver, such as base transceiver 40. 

shown at reference numeral 300. Voltage samples 300 must 5 This process of calculating a spatially averaged radio 

be taken at a rate greater than the bandwidth of the signal frequency signature for implicitly measured RF signal 

being measured. propagation characteristics may be more completely under- 

To calculate impulse response vector h(i), groups of stood with reference to FIG. 7. 

voltage samples s(i), such as sample groups 302 and 304 In FIG. 7, the process of calculating a spatially averaged 

having elements v-v^.j, are input into impulse response 10 radio frequency signature for the selected signature region 

estimator 306. Impulse response estimator 306 delivers from implicitly measured radio frequency signal propaga- 

samples of the estimated channel impulse response, as in tion characteristics begins at block 700 and thereafter passes 

equation 3, and h(i) in FIG. 3. The channel impulse response to block 702, wherein the process selects a first location 

may be estimated according to the method described in within the selected signature region where radio frequency 

"Joint Angle and Delay Estimation Using Shift-Invariance 15 signal propagation characteristics were measured. 

Techniques," IEEE Transactions on Signal Processing, Feb- Next, the process recalls the measured channel impulse 

ruary 1998. responses for the selected location, as illustrated at block 

In a preferred embodiment, when multiple antennas are 704. These channel responses are preferably represented as 

used, an impulse response is measured for each antenna, vectors that are associated with the location at which they 

wherein the impulse response for antenna e is denoted by 20 were collected. If more than one antenna is measured, the 

h c (i). For example, if E number of antennas are used, the impulse responses are recalled as a matrix with each column 

impulse responses for all antennas for each measurement representing the impulse response from a single antenna, as 

may be arranged in matrix H(i) as shown in equation 4. shown in equation 4. 

Next, the process vectorizes the recalled impulse 

lHi>[h.(i) h,(n . . . HJi)] equation 4 , r . . 1( . t . r . 

25 responses and appends the resulting vector as a column in 

After the radio frequency signal propagation characteris- matrix F, wherein the column corresponds to the selected 

tics are measured, the process determines whether or not location, as depicted at block 706. Note that the vectorizing 

there is a next location to measure in the communication process is only necessary if more than one antenna is used 

system coverage area, as depicted at decision block 206. If at the base transceiver to measure channel impulse 

there is another location to measure, the process passes to 30 responses. For L measurement locations, the matrix F can be 

block 208, wherein the next location is selected. After represented as follows: 

selecting the next location, new measurements are taken, as j{i)**vcc{ii(0} 
described above with reference to block 204. Note that the 

selected locations need not be unique, which might be the H^ 1 )/! 2 ) • • 

case if one drives or walks through an intersection in the 35 Next, the process determines whether or not all the 

coverage area twice. Even though more than one impulse measured impulse responses for the measured locations 

response vector may be stored for the same selected within the selected signature region have been appended in 

location, several different locations should be selected matrix F, as illustrated at block 708. If the responses from all 

within a signature region. the locations have not been recalled and vectorized, the 

If, as determined at block 206, all locations have been 40 process selects the next location, as depicted at block 710, 

measured, the process then defines a plurality of signature and returns to block 704 to recall the next measured channel 

regions, wherein each signature region contains a plurality impulse response. 

of selected locations, as illustrated at block 210. In a If the process determines at block 708 that all data has 

preferred embodiment, the selected signature regions are been recalled and appended to form matrix F, the process 

adjacent to each other, and are geometrically similar in shape 45 continues at block 712, wherein principle subspace T is 

and size, so that the regions look like a grid of tiles laid over extracted from matrix F. This principle subspace extraction 

the communications system coverage area. As an example, is generally accomplished by partitioning the result of either 

signature region 44 in FIG. 1 contains a plurality of selected a singular value decomposition operating directly on matrix 

locations 42 from which RF signal propagation character- F, or by an eigendecomposition operating on a sample 

istics have been measured. 50 covariance matrix estimated from data matrix F. The parti - 

In an alternative embodiment, signature regions 44 may tioning separates significant terms from less important terms 

be selected based upon local correlations observed in the that may be heavily distorted by noise or other artifacts. This 

channel impulse responses as determined by comparing a partitioning may be simply accomplished by ranking the 

measurement atone location with a measurement at another, computed singular values or eigenvalues in order of their 

or by comparing a group of measurements with another 55 magnitude, and selecting the P^N largest values from the 

group of measurements. By performing this analysis in the group of N values, or components may be selected that 

defining of signature regions, irregularly shaped signature exceed a threshold value. More sophisticated methods of 

regions may be formed, and those signature regions may determining the principle subspace rank are well known in 

more accurately represent radio frequency propagation char- the literature, e.g. as described in Simon Haykin, ed., 

acteristics from the region, which in turn will make location 60 Advances in Spectrum Analysis and Array Processing, Vol. 

estimation more accurate. ///, Prentice Hall, 1995, pp.5-9 and pp. 21-23. 

Once the signature regions are defined, the process selects After the step represented at block 712 has been 

a first signature region, as depicted at block 212. For the completed, matrix T becomes the spatially averaged RF 

selected signature region, each of the measured RF signal signature for the selected signature region. Thus, the process 

propagation characteristics is processed to produce a spa- 65 of producing a radio frequency signature for the selected 

tially averaged radio frequency signature for the signature signature region from implicitly measured radio frequency 

region, as illustrated at block 214. This spatially averaged signal propagation characteristics ends at block 714. 
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The process described in relation to FIG. 7 is one that 
extracts a principle subspace from implicitly measured radio 
frequency propagation characteristics. For a more complete 
discussion of principle subspace extraction on implicitly 
measured radio frequency propagation characteristics, along 5 
with other methods of extracting principle subspace, see the 
dissertation entitled: "Estimation of Parametric Channel 
Models in Wireless Communication Networks", by 
Michaela C. Vanderveen, May 1998, Stanford University. 

The principal subspace extraction captures the dominant ao 
features of the RF signal propagation characteristics from 
the signature region. This process filters the information 
received and recorded implicitly in the plurality of impulse 
response vectors associated with the selected signature 
region and removes information deemed to be unimportant is 
or less important to location finding. 

An alternative method of extracting the dominant 
features, and one that is conceptually easier to understand, is 
to nonlinearly average the measured impulse responses 
associated with a signature region to deliver a single average 20 
power delay profile. With reference to FIG. 10, the measured 
impulse responses can be converted into power delay pro- 
files by taking the magnitude of the complex valued samples 
of the impulse response at each delay. Plot 1000 is a 
graphical representation of a collection of such power delay 25 
profiles, with power represented in decibels, delay in 
microseconds, and the location of each measurement loosely 
represented as an integer index. The collection of profiles in 
plot 1000 advance along location axis 1002, wherein the 
ofikets along location axis 1002 represent the plurality of 30 
measurements collected that lie within a selected signature 
region. Averaging along location axis 1002 allows us to 
extract the mean features, yielding an average power delay 
profile. With reference to FIG. 11, plot 1100 is a graphical 
representation of such an average power delay profile 35 
extracted from the data of FIG. ,10. The average power delay 
profile can be computed per each antenna measured. The 
average power delay profile, like the principle subspace, can 
be used to represent the average or dominant features of the 
radio frequency signal propagation characteristics. 40 

The dominant features of the radio frequency signal 
propagation characteristics will tend to be associated with 
the location of the transmitter and receiver. Because the base 
transceiver is typically stationary in a traditional wireless 
communications system, we can use the dominant features 45 
to identify the location of the subscriber unit. 

After the spatially averaged radio frequency signature is 
produced for a selected signature region, the process records 
an association between the spatially averaged radio fre- 
quency signature and a location representing the selected 50 
signature region, as depicted at block 216. Preferably this 
recording is done by writing a record in a database that 
associates each spatially averaged radio frequency signature 
with its respective signature region, as represented by a 
location, such as location 46 (See FIG. 1) at the center of 55 
signature region 44. 

Next, the process determines whether or not another 
spatially averaged radio frequency signature needs to be 
, calculated for another signature region in the coverage area, 
as illustrated at block 218. If additional processing for 60 
another signature region is needed, the process selects the 
next signature region, as illustrated at block 220, and then 
returns to block 214 wherein the next spatially averaged 
radio frequency signature is calculated. 

If the process determines at block 218 that all the signa- 65 
lure regions have been processed, the process of database 
formation ends, as depicted at block 222. The database 



,506 Bl 

8 

contains spatially averaged radio frequency signatures for 
each signature region. The spatially averaged radio fre- 
quency signatures contain information that represents the 
dominant radio frequency propagation characteristics for 
signals transmitted from the signature region, as determined 
by mathematically filtering less important information from 
the measurements from a plurality of selected locations in 
the signature region. The database is now ready to use to 
estimate a subscriber unit location, as described in detail 
below. 

With reference now to FIG. 4, there is depicted a high- 
level logic flowchart that illustrates the method and opera- 
tion of measuring a time averaged radio frequency signature 
for a transmitter in a coverage area of a wireless commu- 
nication system, comparing the time averaged radio fre- 
quency signature with a spatially averaged radio frequency 
signature, and in response to the comparison, generating a 
figure of merit in accordance with the method and system of 
the present invention. A collection of figures of merit can 
then be used to estimate the transmitter's location by decid- 
ing from which signature region the transmitter is likely to 
be transmitting. 

As illustrated, the process begins at block 400, and 
thereafter passes to block 402 wherein the process receives 
a request to estimate the location of a subscriber unit. 

After receiving the request, the process measures the radio 
frequency signal propagation characteristics of a signal from 
the transmitter of the subscriber unit for which the system is 
estimating a location, as illustrated at block 404. The mea- 
sured radio frequency signal propagation characteristics 
may include the number of multipath components, the 
average power in each multipath component, the gain and 
phase of the antenna due to the azimuth and elevation angles 
of arrival and polarization of the wave front, the Doppler 
frequency of the wave front, and the relative time delay of 
the wave front. These characteristics may be measured 
explicitly by parametric decomposition of the received sig- 
nal. However, in a preferred embodiment, the radio fre- 
quency signal propagation characteristics are measured 
implicitly by sampling and recording the channel impulse 
responses of the particular channel between the subscriber to 
be located and a base transceiver. 

The measurement of the radio frequency signal propaga- 
tion characteristics of the signal from the subscriber unit to 
be located is very similar to measuring the radio frequency 
signal propagation characteristics of a signal from a trans- 
mitter at a selected location during the formation of the 
spatially averaged radio frequency signature database, as 
illustrated at block 204 in FIG. 2. 

Once radio frequency signal propagation characteristics 
have been measured, the process determines whether or not 
another measurement is needed, as depicted at block 406. 
Preferably, several measurements of RF signal propagation 
characteristics will be taken periodically so that a time 
averaged radio frequency signature for the subscriber unit to 
be located can be calculated. A single instantaneous mea- 
surement of the RF signal propagation characteristics may 
be adequate in some circumstances, but it is well known that 
a single sample estimator can be subject to large variations 
due to the underlying random processes in signal propaga- 
tion. Estimating the RF signal propagation characteristics 
from multiple measurements uses statistical principles to 
mitigate the randomness and extract the average structure or 
dominant features of the RF signal propagation character- 
istics. 

If another measurement is needed, the process delays, as 
depicted at block 408, and then passes back to block 404 
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wherein radio frequency signal propagation characteristics After the step represented at block 512 has been 

are measured again. The number of measurements taken completed, matrix S becomes the time averaged RF signa- 

may be determined by the number of periodic measurements ture for the signal from the subscriber unit to be located, 

that can be taken over a predetermined fixed interval. Thus, the process of producing a time averaged radio 

If the process determines at block 406 that enough mea- 5 frequency signature from implicitly measured radio fre- 

surements have been taken, the process proceeds to block quency signal propagation characteristics ends at block 514. 

410 wherein the measured radio frequency signal propaga- The process described in relation to FIG. 5 is one that 

tion characteristics are processed to produce a time averaged extracts principle subspace from implicitly measured radio 

radio frequency signature for the signal from the subscriber frequency propagation characteristics. For more complete 

unit to be located, wherein the time averaged radio fre- 10 discussion of principle subspace extraction on implicitly 

quency signature represents the dominant features of the measured radio frequency propagation characteristics, along 

radio frequency signal propagation characteristics. This pro- with other methods of extracting principle subspace, see the 

cess of calculating the time averaged radio frequency sig- dissertation entitled: "Estimation of Parametric Channel 

nature may be more completely understood with reference to Models in Wireless Communication Networks", by 

FIG. 5. is Michaela C. Vanderveen, May 1998, Stanford University. 

In FIG. 5, the process of calculating a time averaged radio The principal subspace extraction captures the dominant 

frequency signature from implicitly measured radio fre- features of the signal propagated from the signature region, 

quency signal propagation characteristics begins at block This process filters the information received and recorded 

500 and thereafter passes to block 502 wherein the process implicitly in the plurality of impulse response vectors asso- 

selects a first time at which radio frequency signal propa- 20 ciated with the selected signature region and removes infor- 

gation characteristics were measured for the signal from the mation deemed to be unimportant or less important to 

subscriber unit to be located. location finding. 

Next, the process recalls the measured channel impulse An alternative method of extracting the dominant 
responses for the selected time, as illustrated at block 504. features, and one that is conceptually easier to understand, is 
These channel impulse responses are preferably represented 25 to nonlinearly average the measured impulse responses 
as vectors that are associated with the time at which they associated with a signature region to deliver a single average 
were collected. If more than one antenna is measured, the power delay profile. With reference to FIG. 8, the measured 
responses are recalled as a matrix with each column repre- impulse responses can be converted into power delay pre- 
senting the impulse response from a, single antenna, as in files by taking the magnitude of the complex valued samples 
equation 4. 30 of the impulse response at each delay. Plot 800 is a graphical 

Next, the process vectorizes the recalled impulse representation of a collection of such power delay profiles, 

responses and appends the resulting vector as a column in with power represented in decibels, delay in microseconds, 

matrix G, wherein the column corresponds to the selected and the time of measurement of each power delay profile 

time, as depicted at block 506. Note that the vectorizing measured in milliseconds. The collection of profiles in plot 

process is only necessary if more than one antenna is used 35 800 advances along time axis 802, wherein location axis 802 

at the base transceiver to measure channel impulse represents the plurality of measurements collected within a 

responses. For L measurement times, the matrix G can be span of time. Averaging along time axis 802 allows us to 

represented as follows: extract the mean features, yielding an average power delay 

oC\m (hc\\ profile. With reference to FIG. 9, plot 900 is a graphical 

40 representation of such an average power delay profile 

<M*9i) s@) ■ • • s(L)] extracted from the data of FIG. 8, The average power delay 

Next, the process determines whether or not all the profile can be computed per each antenna measured. The 

measured impulse responses for all the measured times have average power delay profile, like the principle subspace, can 

been appended in matrix G, as illustrated at block 508. If all be used to represent the average or dominant features of the 

the times have not been recalled, vectorized, and appended 45 radio frequency signal propagation characteristics, 

into G, the process selects the next time, as depicted at block The dominant features of the radio frequency signal 

510, and returns to block 504 to recall the next measured propagation characteristics will tend to be associated with 

channel impulse response. the location of the transmitter and receiver. Because the base 

If the process determines at block 508 that all data has transceiver is typically stationary in a traditional wireless 

been recalled and appended to form matrix G, the process 50 communications system, we can use the dominant features 

continues at block 512, wherein principle subspace S is to identify the location of the subscriber unit, 

extracted from matrix G. This principle subspace extraction With reference again to FIG, 4, after the time average 

is generally accomplished by partitioning the result of either radio frequency signature has been calculated at block 410, 

a singular value decomposition operating directly on matrix the process selects a data record from the signature database 

G, or by an eigendecomposition operating on a sample 55 having data that represents the dominant features of the 

covariance matrix estimated from data matrix G. The par- radio frequency signal propagation characteristics of a 

titioning separates significant terms from less important selected signature region, as illustrated at block 412. 

terms, which may be heavily distorted by noise or other The process then compares the time averaged radio fre- 

artifacts. This partitioning may be accomplished by ranking quency signature for the signal from the subscriber unit to be 

the computed singular values or eigenvalues in order of their 60 located to the selected data record from the signature 

magnitude, and selecting the P^N largest values from the database, which contains a spatially averaged radio fre- 

group of N values. Alternatively, components may be quency signature, to produce a figure of merit, as depicted at 

selected that exceed a threshold value. More sophisticated block 414. In a preferred embodiment, where the radio 

methods of determining the principle subspace rank are well frequency signatures are derived from implicit 

known in the literature, e.g. as described in Simon Haykin, 65 measurements, the process of computing a figure of merit 

ed., Advances in Spectrum Analysis and Array Processing, may be accomplished by calculating the distance between 

Vol III, Prentice Hall, 1995, pp.5-9 and pp. 21-23. the principle subspaces represented by the time averaged 
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radio frequency signature from the subscriber unit to be Doppler frequency of the wave front, and the relative time 

located and the spatially averaged radio frequency signature delay of the wave front. As mentioned previously, these 

from the selected data record in the database. The distance characteristics may measured explicitly by parametric 

computation may be performed by computing the norm of decomposition of the received signal, or implicitly by 

the projection of one signature onto the other. This com- 5 estimating, sampling and recording the channel impulse 

parison may be a weighted comparison where more impor- responses. 

tant information is given more consideration. The output of channel estimator 612, whether implicit 

Alternatively, for explicitly measured radio frequency representations of radio frequency propagation characteris- 

signal propagation characteristics, the signatures may be tics or explicit measurements of these characteristics, are 

hyperdimensional vectors representing the parameters of the 10 stored in storage 614. This storage may be implemented with 

explicit propagation model, as in equation 2, and the figure any commonly known data storage device, such as a hard 

of merit representing their comparison can be the distance disk, tape, random access memory, CD ROM, and the like, 

between spatially and temporally averaged versions of these Processor 616 is coupled to storage 614 for recalling and 

hyperdimensional vectors, similar to the methods outlined processing data, and storing the results of such processing, 

above for implicit measurements. Note that here, as 15 Processor 616 may be implemented with a general purpose 

elsewhere, averaging can include traditional linear microprocessor, a custom application specific integrated 

averaging, nonlinear averaging, or subspace based averag- circuit, or a digital signal processor. Within processor 616 

ing via extraction of a principle subspace. are other functional blocks which may be implemented with 

Next, the process stores the figure of merit associated with various combinations of hardware and software. Note also 

the selected data record from the spatially averaged radio 20 that channel estimator 612 may also be implemented in 

frequency signature database, as illustrated at block 416. processor 616. 

Thereafter, the process determines whether or not there is With regard to a system that builds the spatially averaged 

another data record to compare, as depicted at block 418. If radio frequency signature database, processor 616 contains 

there is another data record to compare, the process selects region selector 618, which recalls data stored in storage 614, 

the next data record, as illustrated at block 420, and then 25 wherein the data represents measured radio frequency signal 

returns to block 414 to compare the two radio frequency propagation characteristics of a particular signature region, 

signatures. Region selector 618 then passes this data to signature region 

If the process determines that all data records have been processor 620, which computes a spatially averaged radio 

compared, the process continues at block 422, wherein the frequency signature for a selected signature region and 

process ranks the data records based on the associated 30 stores such signature in database 622, where it is associated 

figures of merit in order to find a location estimate. Once the with a particular signature region. The spatially averaged 

location estimate has been found, the process outputs the radio frequency signatures are calculated as discussed with 

coordinates of a signature region that corresponds to the reference to FIG. 2 and 7. 

location estimate, as depicted at block 424. This output may When processor 616 has received a request to estimate the 

be the coordinates of a location associated with the signature 35 location of a subscriber unit, subscriber signature processor 

region, such as the coordinates of location 46. The process 624 recalls data measurements made by channel estimator 

then terminates as indicated at block 426 . 612 from storage 614. These measurements represent radio 

With reference now to FIG. 6, there is depicted a system frequency signature propagation characteristics of a signal 

used to create a radio frequency signature for a signature from the subscriber unit to be located measured over a 

region in a coverage area of a wireless communication 40 plurality of times. Subscriber signature processor 624 com- 

system, and also for measuring a time averaged radio putes a time averaged radio frequency signature, as 

frequency signature for a subscriber unit to be located in the described with reference to FIGS. 4 and 5. 

coverage area, in accordance with a method and system of The process of comparing a time averaged radio fre- 

the present invention. As illustrated, a plurality of wave quency signature with a spatially averaged radio frequency 

fronts 600 and 602 are received by antennas 604-606. 45 signature begins with signature region searcher 626 recalling 

Wavefront 600-602 represent a signal from a single source a spatially averaged radio frequency signature from database 

that has been subjected to multipath distortion. In other 622 for a selected signature region. The output of signature 

words, both wave front 600 and 602 came from the same region searcher 626 is coupled to one of the inputs of 

transmitter through different propagation paths. signature comparator 628, and another input of signature 

In a typical embodiment, antennas 606-604 may be 50 comparator 628 is coupled to the output of subscriber 
implemented with two diversity receive antennas in each of signature processor 624. Signature comparator 628 corn- 
three sectors on the top of a typical cellular tower. Antennas pares the time averaged radio frequency signature of the 
606-604 are coupled to receiver and down converter 608, signal from the subscriber unit to be located and the spatially 
which translates the spectrum from radio frequency to averaged radio frequency signature of a selected signature 
complex valued baseband. 55 region and generates a figure of merit as a result of the 

The output of receiver and down converter 608 is coupled comparison. Signature comparator 628 may be coupled to 

to the input of analog- to-digital converter 610, which con- storage 614 in order to temporarily store figures of merit that 

verts the complex valued baseband signal to a stream of relate to each signature region recalled by signature region 

digital values. The stream of digital values output by analog- search 626. Signature comparator 628 may also store the 

to-digital converter 610 is coupled to the input of channel 60 figure of merit and the corresponding signature region in a 

estimator 612. Channel estimator 612 estimates the radio peak holding manner such that if a higher figure of merit is 

frequency signal propagation characteristics of wave fronts produced, that figure of merit, along with its signature region 

600-602. These radio frequency signal propagation charac- coordinates, replace the previous value. Signature compara- 

teristics may include the number of multipath components, tor 628 is coupled to output function 630, which formats the 

the average power in each multipath component, the gain 65 subscriber unit's estimated location represented by the sig- 

and phase of the antenna due to the azimuth and elevation nature region with the best figure of merit. This output may 

angles of arrival and polarization of the wave front, the be in the form of position coordinates of the signature 



03/10/2004, EAST Version: 1.4.1 



US 6,256,506 Bl 

13 14 

region, such as location 46, or by displaying a map showing spatially averaged radio frequency signature for the selected 

the location of the subscriber unit. signature region further includes decomposing a structured 

It should be noted that the above embodiments have matrix containing vector representations of the signal propa- 

tended to speak to an uplink mode of operation, in which a gation characteristics within the selected signature region to 

mobile subscriber unit is the transmitter, and the base 5 extract a subspace that represents the dominant features of 

transceiver station is the receiver. Alternate combinations of the rad i 0 frequency signal propagation characteristics in the 

receiver and transmitter are possible in practice, including selected signature region. 

operation on the downlink. Downlink operation may be 4 Jhe method for creating a radio frequency signature 

achieved by measuring the radio frequency signal propaga- according t0 claim 3 wherein the step of decomposing a 

tion characteristics in the subscriber unit, and either com- 1Q structured matrix cont aining vector representations of the 

puting an associated radio frequency signature in the sub- & . { ropagation characteristics further includes the steps 

scriber unit, or sending the characteristics to the base 0 £. 

transceiver for processing to a radio frequency signature. * 

Similar to the uplink methods, this can be done for multiple performing an eigendecomposition on the structured 

antennas at the base transceiver. „ matnx staining vector representations of the sign a 

The foregoing description of a preferred embodiment of 15 Propagation characteristics to produce a plurality of 

the invention has been presented for the purpose of illus- eigenvector and their corresponding eigenvalues; and 

tration and description. It is not intended to be exhaustive or selecting N of the eigenvectors that correspond to the N 

to limit the invention to the precise form disclosed. Obvious largest eigenvalues. 

modifications or variations are possible in light of the above 5. The method for creating a radio frequency signature 
teachings. The embodiment was chosen and described to according to claim 1 wherein the step of measuring radio 
provide the best illustration of the principles of the invention frequency signal propagation characteristics associated with 
and its practical application, and to enable one of ordinary a signal transmitted from a plurality of locations in the 
skill in the art to utilize the invention in various embodi- coverage area of the wireless communication system further 
ments and with various modifications as are suited to the includes measuring radio frequency signal propagation char- 
particular use contemplated. All such modifications and acteristics in the form of an impulse response for a plurality 
variations are within the scope of the invention as deter- of locations in the coverage area of the wireless communi- 
mined by the appended claims when interpreted in accor- cation system. 

dance with the breadth to which they are fairly, legally, and * The method for creating a radio frequency signature 

equitably entitled. according to claim 5 wherein the step of implicitly measur- 

What is claimed is: 30 m S radio frequency signal propagation characteristics of an 

1. A method for creating a radio frequency signature for impulse response for a plurality of locations in the coverage 
a signature region in a coverage area of a wireless commu- area further includes measuring a power delay profile for a 
nication system, the method comprising the steps of: plurality of locations in the coverage area of the wireless 

measuring radio frequency signal propagation character- 3S communication system. 

istics associated with a signal transmitted from a plu- 7 - A s y slem for creatin g a rad, ° frequency signature for a 

rality of locations in the coverage area of the wireless signature region in a coverage area of a wireless communi- 

communication system; cation svstem comprising: 

defining a plurality of signature regions in the coverage means for measuring radio frequency signal propagation 

area that contain selected locations having selected 40 characteristics associated with a signal transmitted 

measured radio frequency signal propagation charac- frora a plurality of locations in the coverage area of the 

teristics* wireless communication system; 

for a selected signature region, processing the selected means for defining a plurality of signature regions in the 

measured radio frequency signal propagation charac- coverage area that contain selected locations having 

teristics within the selected signature region to produce 45 selected measured radio frequency signal propagation 

a spatially averaged radio frequency signature for the characteristics; 

selected signature region, wherein the processing means for processing the selected measured radio fre- 
extracts the dominant features of the radio frequency quency signal propagation characteristics within a 
signal propagation characteristics of a signal transmit- selected signature region to produce a spatially aver- 
ted within the selected signature region; and 50 aged radio frequency signature for the selected signa- 
recording an association between the spatially averaged [ ure region, wherein the processing extracts the domi- 
radio frequency signature and the selected signature riant features of the radio frequency signal propagation 
region, characteristics of a signal transmitted within the 

2. The method for creating a radio frequency signature selected signature region; and 

according to claim 1 wherein the step of defining a plurality 55 means for recording an association between the spatially 

of signature regions in the coverage area that contain averaged radio frequency signature and the selected 

selected locations having selected measured radio frequency signature region. 

signal propagation characteristics further includes defining a 8. The system for creating a radio frequency signature 

plurality of signature regions in the coverage area that according to claim 7 wherein the means for defining a 

contain selected locations having selected measured radio 60 plurality of signature regions in the coverage area that 

frequency signal propagation characteristics, wherein the contain selected locations having selected measured radio 

signaiure regions are adjacent to each other, and are geo- frequency signal propagation characteristics further includes 

metrically similar in shape and size. means for defining a plurality of signature regions in the 

3. The method for creating a radio frequency signature coverage area that contain selected locations having selected 
according to claim 1 wherein the step of processing the 65 measured radio frequency signal propagation 
selected measured radio frequency signal propagation char- characteristics, wherein the signature regions are adjacent to 
acteristics within the selected signature region to produce a each other, and are geometrically similar in shape and size. 
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9. The system for creating a radio frequency signature 
according to claim 7 wherein the means for processing the 
selected measured radio frequency signal propagation char- 
acteristics within the selected signature region to produce a 
spatially averaged radio frequency signature for the selected 5 
signature region further includes means for decomposing a 
structured matrix containing vector representations of the 
signal propagation characteristics within the selected signa- . 
ture region to extract a subspace that represents the domi- 
nant features of the radio frequency signal propagation 10 
characteristics in the selected signature region. 

10. The system for creating a radio frequency signature 
according to claim 9 wherein the means for decomposing a 
structured matrix containing vector representations of the 
signal propagation characteristics further includes: 15 

means for performing an eigendecomposition on the 
structured matrix containing vector representations of 
the signal propagation characteristics to produce a 
plurality of eigenvectors and their corresponding eigen- 
values; and 
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means for selecting N of the eigenvectors that correspond 
to the N largest eigenvalues. 

11. The system for creating a radio frequency signature 
according to claim 7 wherein the means for measuring radio 
frequency signal propagation characteristics associated with 
a signal transmitted from a plurality of locations in the 
coverage area of the wireless communication system further 
includes means for measuring radio frequency signal propa- 
gation characteristics in the form of an impulse response for 
a plurality of locations in the coverage area of the wireless 
communication system. 

12. The system for creating a radio frequency signature 
according to claim 11 wherein the means for implicitly 
measuring radio frequency signal propagation characteris- 
tics of an impulse response for a plurality of locations in the 
coverage area further includes means for measuring a power 
delay profile for a plurality of locations in the coverage area 
of the wireless communication system. 
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